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First results of a detailed analysis of p+p elastic scattering data are 
presented from ISR to LHC energies utilizing the quark-diquark model of 
protons in a form proposed by Bialas and Bzdak. The differential cross- 
section of elastic proton-proton collisions is analyzed in detailed and sys- 
tematic manner at small momentum transfers, starting from the energy 
range of CERN ISR at = 23.5 GeV, including also recent TOTEM 
data at the present LHC energies at y/s = 7 TeV. These studies confirm 
the picture that the size of protons increases systematically with increas- 
ing energies, while the size of the constituent quarks and diquarks remains 
approximately independent of (or only increases only slightly with) the col- 
liding energy. The detailed analysis indicates correlations between model 
parameters and also indicates an increasing role of shadowing at LHC en- 
ergies. 
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1. Introduction 

The differential cross section of elastic scattering of p+p collisions allows 
one to study the internal structure of protons using the theory of diffraction. 
Varying the momentum transfer one can change the resolution of the inves- 
tigation: increasing the momentum transfer corresponds to looking more 
and more deeply inside the structure of protons. One of the fundamental 
outcomes of diffractive p+p scattering studies was the indication that pro- 
tons have a finite size and a complicated internal structure, thus the protons 
can be considered as composite objects. 

Our interest in this problem has been triggered by two factors: an inter- 
esting series of recent theoretical work and also new data from the TOTEM 
experiment at CERN LHC. These are detailed below. Recently, we became 
aware of an inspiring series of papers of Bialas, Bzdak and collaborators, 
who studied elastic proton-proton p[], pion-proton ^ and nucleus- nucleus 
collisions t3],[l] in a framework where the proton was considered as a com- 
posite object that contains correlated quark and diquark constituents. In 
this work, we confirm their main conclusion: the quark-diquark model of nu- 
cleon structure at low momentum transfer does capture the main features 
of this problem and indeed it deserves a closer, more detailed attention. 
This study is dedicated to a follow-up, more detailed investigation, not only 
including an estimation of the best values of the parameters, as was done 
in ref. [1], but also determining their errors and also the evaluation of the 
models based on an analysis of the presented fit quality. In order to reach 
these goals, we utilized standard experimental techniques, such as multi- 
parameter optimalization or fitting with the help of the MINUIT function 
minimalization and multi-parameter optimalization package [5]. 

In addition to these simple, straightforward and interesting theoretical 
investigations of elastic scattering data from CERN ISR in the energy range 
of ^/s = 23.5, 30.7, 52.9 and 62.5 GeV, that were already analyzed in ref. [T], 
new elastic scattering data became available recently at ^/s = 7 TeV [TB] 
from the CERN LHC experiment TOTEM. So we have tested the model 
of Bialas and Bzdak not only at ISR energies but also at the currently 
available highest LHC energies on recent TOTEM data, in order to learn 
more details about the evolution of the properties of p+p elastic interactions 
in the recently opened, few TeV energy range. 

Most of the arguments for the composite structure of hadrons have been 
derived from the studies of lepton-hadron interactions. The emerging stan- 
dard picture is that hadrons are either mesons, composed of valence quarks 
and anti-quarks, or (anti)baryons, composed of three valence (anti)quarks, 
that carry the quantum numbers, while the electrically neutral gluons carry 
color charges and provide the binding among the quarks and anti-quarks. 
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The exact contribution of quarks, gluons, and the sea of virtual quark- 
antiquark pairs and gluons to certain hadronic properties e.g. spin is still 
under detailed investigation. For example, the gluon contribution to the pro- 
ton spin is still not fully constrained [6]. Also, more than 10 exotic hadronic 
resonances called X, Y and Z states were recently discovered in electron- 
positron collisions at the world's highest luminosities in the BELLE experi- 
ment at KEK. These hadronic states cannot be interpreted in the standard 
picture of quark-antiquark or three (anti)quark bound states, according to 
refs. [71 IS] Thus even nowadays there are still several open questions that 
are related to the compositeness of the hadrons in general. As gluons do 
not interact with directly with leptons, their properties are best explored 
with the help of the strong hadronic interactions. For example, the gluon 
contribution to the proton spin is investigated with the help of polarized 
proton - polarized proton collisions at RHIC [6j. 

In this work, we focus on the effects of internal correlations between the 
quarks inside the protons, examining in detail proton-proton elastic scatter- 
ing at several ISR and also at the currently maximal available LHC energies. 
Let us recall, that a similar analysis involving three independent quarks was 
not able to properly describe the ISR data [9j. In that model quarks were 
considered as "dressed" valence quarks in the sense that they contain the 
gluonic and qq contribution as well, as if the glue would be concentrated 
around pointlike valence quarks. 30 years after the three independent quark 
model of ref. [S] another three-quark model of the protons was proposed in 
ref. [1], that included interesting correlations between two dressed valence 
quarks to form a diquark. This quark-diquark model [Ij, is the basis of our 
current study, in order to examine in details proton-proton elastic scattering 
at ISR and LHC energies. 

The quark-diquark picture of elastic p+p scattering resembles to the 
Glauber optical model |12j in nuclear physics, where a multiple expansion 
is applied. The Glauber model, developed originally for nuclear multiple 
scattering problems like cross sections of protons and neutrons on deuteron, 
became a standard model of high energy interactions in nuclear physics 
where multiple interactions are built up from superpositions of nucleon- 
nucleon scattering. This model became a fundamental and successfull tool 
in describing nuclear collisions at high energy |10j . 

The body of this manuscript is organized as follows: the theoretical 
models are described in Section [51 for both cases, when the diquark acts as 
a single entity and also when it acts as a composite object. In Section [3] the 
original results of Bialas and Bzdak are reproduced and detailed, including 
a note on the quality of these data descriptions. Section H] contains our 
main new results, the first MINUIT results. The conclusion takes place in 
Section [5l 
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2. Elastic scattering in the quark- diquark model 

We describe proton-proton interactions as collision of two systems, each 
one composed of a dressed quark and diquark. The collision is schematically 
illustrated on Fig. [TJ 

V\ Proton 1 




Fig. 1. The schematic scattering situation of the two protons when the diquark is 
assumed to be a single entity. 

The interaction between quarks and diquarks is assumed to be purely 
absorptive. Consequently the amplitude has no real part and the imaginary 
part - dominating at high energy - is given by the absorption of the incoming 
particle wave, namely the inelastic (non-diffractive) collisions. 

In the impact parameter space the inelastic proton-proton cross-section 
for a fixed impact parameter b can be given by the following formula [T] 

+ 00 +00 

a(b) = j ... j d^Sqd^s'qd^Sd<i^s'^D{sq,Sd)D{sq ,Sd)a{sq,s^^^^ 

— oo — oo 

(1) 

where s^, Sq and s^, s^' are the transverse positions of the quarks and 
diquarks respectively, and the integrand is the probability function of having 
inelastic interaction at a given impact parameter h and transverse positions 
of the constituents. 

The quark-diquark distribution inside the nucleon is taken into account 
with the following Gaussian 

D {s-q, s-d) = ^-±^e-^'l+'l^"'l5\s-'d + Xsq), A = mq/rud, (2) 

where Rp is the "proton size" , the variance of the distribution and A is the 
mass ratio of the quark and the diquark. Obviously 1/2 < A < 1, where 
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1/2 would indicate a loosely bound diquark. The two dimensional delta 
function preserves the center of mass in the tranverse plane. 

Elastic interactions of the constituents are independent inside the proton, 
accordingly the probability distribution of elastic proton-proton collision 
is the product of the probability distribution of elastic interactions of the 
constituents [HI [13] 

C7(s^,Sd;S^',Sd';6) = 1 - Y\ {1 - '7afe(& + Sa' - 4')}- (3) 

a,be{q,d} 

The inelastic differential cross-sections are parametrized with Gaussian dis- 
tributions 

^ab (s1 = Aabe-''/^l», Rl, = Rl + Rl (4) 

where Rab is the variance of having an inelastic collision, which is calculated 
from the sum of the squared Rq, Rd radius parameters; the Aab parameters 
are the ampitudes. From unitarity the elastic amplitude in impact param- 
eter spacqj 



tel{h) = l-^l-<j{b). (5) 

The elastic amplitude in momentum transfer representation is the Fourier- 
transform of the amplitude in impact parameter space 

+00 +00 +00 

r(A)= j j t,i(b)e'^-^d% = 27r J t^i (b) Jq (Ab) bdb, (6) 

— oo — oo 

where A = |A|, b = \b\ and Jq is the zeroth Bessel-function of the first kind. 
Then the elastic differential cross section reads as 

2.1. The diquark is assumed to act as a single entity 

The subject of this section is to analyse that case when the quark and 
diquark radii are independent model parameters, which means that the 
diquark is considered as one entity as indicated on Fig. [TJ In this case, the 
number of free parameters can be reduced if we assume that the number 
of partons is twice as many in the diquark than in the quark. Prom the 
inelastic differential cross sections ^ the total inelastic cross sections are 

+ 00 +00 

crab= j j CTab (sjd^s = irAabRlb, a,b £ {q, d}. (8) 



^ As it was mentioned the real part of the amplitude is ignored. 
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Our assumption tells us that 



c^qq ■ CTqd : (T dd =1:2:4, 



(9) 



from which we can deduce the following expressions 



ml 



Add 



A, 



ml 



(10) 





which means that every A^h parameter can be expressed in term of Aqq. 
With these ingredients the calculation of ([T|) reduces to Gaussian integra- 
tions. The general or master formula for these Gaussian integrals is given 
in the Appendix. 

In the next two subsections the elastic proton-proton data analysis are 
presented in the case when the diquark is assumed to have no internal, 
more detailed structure in elastic collisions. First, we demonstrate with 
plots that we reproduce the results of the original paper at ISR energies 
[1]. This forms a solid basis for imrovement as presented in the subsequent 
parts of our current study. Our MINUIT fits are then presented, utilizing 
same ISR data and then the new results are presented at 7 TeV utilizing 
new data of the TOTEM experiment [l6j . 




2.2. Analysis of the case when the diquark acts as composite object 

The scattering when the diquark is a composite object is illustrated on 
Fig. [2j The quark distribution inside the diquark is supposed to have a 
Gaussian shape 



where Sqi and Sg2 are the transverse quark positions inside the diquark, and 






d'^ - Rd- R-l 



(12) 



is the variance of the quark distribution, calculated from the diquark and 
quark radius parameters. 



Fig. 2. The scattering situation of the two protons when the diquark is assumed 
to be composed of two quarks, and the proton symbohcahy can be written as 

p=(q,(qq))- 



If the diquark has the internal structure (fTT]l then the aqd, (Jdq and add 
inelastic differential cross sections dH) can be calculated from Uqq using an 
expansion analogous to expression ([3]). The results for aqd and add are the 
following [1] 

A A 7?2 „2 1 42 r>2 

(^) = #=4- 13^ _ 4g^, ,-.VH? , (13, 



and 



AA\X -^-^ _ _4|^ 

The relevant formula to calculate the inelastic cross section of eq. ([1]) is again 
obtained from the expressions summarized in the Appendix. 
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3. Reproduction and cross-checks of the fits of Bialas and Bzdak 

3.1. The diquark acts as a single entity 

In this section we reproduce the fitting formula of ref. [Ij, based on the 
best values of the model parameters that they have published. These values 
were obtained in ref. [1] from fitting four essential elements of the elastic 
scattering cross-sections: they adjusted the model parameters so that (i) 
the total inelastic cross-section, (ii) the slope of the differential inelastic 
cross-section at t = 0, (iii) the position of the diffractive minimum, and 
(iv) the height of the first diffractive maximum just after the minimum be 
in agreement with the data. This method lead to remarkable simplicity 
in the fitting strategy and a very nice, apparent overall agreement with 
the measured data, as can be seen from the Figures of ref. [1] . As a 
consequence of this method, the errors of the fit parameters in ref. were 
not presented, and the overall fit quality parameters (x^/NDF, CL) were 
not specified as well, although these are the qualifyers that determine the 
acceptability of certain models when the language of mathematical statistics 
is utilized to characterize the data description. By recalculating the original 
curves, we check in this section the quality of our reproduction of the results 
of Bialas and Bzdak, and also supplemented their paper with the previously 
undetermined fit quality parameters. 

The results are shown on Figs. [HHSlfor the energies ^/s =23.5, 30.7, 52.8, 
62.5 GeV respectively. More recently the elastic proton-proton da/dt was 
measured by the TOTEM experiment at LHC. The first TOTEM data were 
published in the |t| range from 0.36 GeV up to 2.5 GeV [l6j. The fit quality 
parameters presented here are restricted to those data points which fall into 
this \t\ region, in order to compare the accuracy of the model in the ISR 
and LHC energy regimes. 
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p+p p+p, diquark as a single entity,Vs=23.5 GeV 



Reproduction of [1] 
Fit range = 0.36 < t < 2.50 GeV^ 
yj- 1 NDF = 674.06/44= 15.32 
CL = 0.000 % 



o Data points 

• Not used 

Theory Ret. [1] 

Extrapolation 




Fig. 3. Reproduction of the calculation of A. Bialas and A. Bzdak from [1 at 23.5 
GeV, in the case when the diquark is assumed to be a single entity. We have 
fixed the values of the model parameters to values given in Ref. [T]. In order to 
determine fit quality, we restricted the fit range to 0.36 GeV < -t < 2.5 GeV, to 
allow a systematic comparision with more recent TOTEM data at LHC. As the 
model of Bialas and Bzdak is known to be singular at the dip, in this and all other 
data analysed we have left out 3 data points from the dip region. 
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p+p p+p, diquark as a single entity,Vs=30.7 GeV 
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Fig. 4. Same as Fig. H but at Vs=30.7 GeV. 





o Data points 
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p+p p+p, diquark as a single entity,Vs=52.8 GeV 



Reproduction of [1] 
Fit range = 0.36 < t < 2.50 GeV^ 
X^/NDF = 105.11/33= 3.19 
CL = 0.000 % 



o Data points 

• Not used 

Theory Ret. [1] 

Extrapolation 




-t [GeV^] 

Fig. 5. Same as Figs. [21&I1 at Vs=52.8 GeV. 
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p+p p+p, diquark as a single entity,Vs=62.5 GeV 
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Reproduction of [1] 
Fit range = 0.36 < t < 2.50 GeV' 
X^/NDF = 102.92/33= 3.12 
CL = 0.000 % 



o Data points 

Not used 
Theory Ret. [1] 
Extrapolation 
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Fig. 6. Same as Figs. except that the energy is -^5=62. 5 GeV. 
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3.2. Reproduction of earlier results for the case of composite diquarks 

We checked the correctness of our implementation of the formula from 
the original paper [1]. In the original analysis the errors of the fit parameters 
were not shown and the overall fit quality parameters (x^/NDF, CL) are 
missing as well. The missing fit quality parameters will be supplemented 
in this section. The results are shown on Fig. [TlfTOl for the energies 23.5, 
30.7, 52.8, 62.5 GeV respectively. More recently the elastic proton-proton 
da/dt was measured by the TOTEM experiment at LHC. The first TOTEM 
data covers the |t| range from 0.36 GeV up to 2.5 GeV [l6j. The fit quality 
parameters presented here are restricted to those data points which fall into 
this region, in order to compare the accuracy of the model in the ISR and 
LHC energy regime. 
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p+p p+p, diquark as a qq entity,\fs=23.5 GeV 



Reproduction of [1] 
Fit range = 0.36 < t < 2.50 GeV^ 
yj- 1 NDF = 1043.64/44= 23.72 
CL = 0.000 % 



Data points 
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Fig. 7. Reproduction of the calculation of A. Bialas and A. Bzdak [T] at 23.5 
GeV, in the case when the diquark is assumed to be a qq entity. We have fixed the 
values of model parameters to values given in [T] . In order to determine fit quality, 
we restricted the fit range to 0.36 GeV < -t < 2.5 GeV, to allow a systematic 
comparision with more recent TOTEM data at LHC. As the model of Bialas and 
Bzdak is known to be singular at the dip, in this and all other data analysed we 
have left out 3 data points from the dip region. 



p+p p+p, diquark as a qq entity,\fs=30.7 GeV 
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Reproduction of [1] 
Fit range = 0.36 < t < 2.50 GeV^ 
X^/NDF = 182.60/32= 5.71 
CL = 0.000 % 
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Same as Fig [71 except that the analyzed energy is ^/s = 30.7 GeV 
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p+p p+p, diquark as a qq entity,\fs=52.8 GeV 



Reproduction of [1] 
Fit range = 0.36 < t < 2.50 GeV^ 
X^/NDF = 331 .11/33= 10.03 
CL = 0.000 % 



o Data points 

Not used 
Theory Ret. [1] 
Extrapolation 
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Fig. 9. Same as Fig El given at = 52.8 GeV. 



p+p p+p, diquark as a qq entity,\fs=62.5 GeV 



Reproduction of [1] 
Fit range = 0.36 < t < 2.50 GeW^ 

I NDF = 206.87/33= 6.27 
CL = 0.000 % 
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Fig. 10. Identical with Fig [71 except that it is evaluated at ^/s = 62.5 GeV 
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4. MINUIT fit results for ISR and TOTEM data 
4.1. The diquark acts as a single entity 

In this section the MINUIT fit results are presented for the ISR |141 
[15] and TOTEM [16j proton-proton elastic scattering data using the single 
entity model to describe the diquark. Results are illustrated on Fig. [TTl 
[T21 The confidence levels, and model parameters together with their errors 
are presented in Table [TJ amended with the calculated total elastic cross 
sections including their uncertainty evaluated from the MINUIT fits. The 
ratios of the inelastic cross sections were fixed with expression ([9]) in order 
to decrease the number of free parameters, therefore these ratios will be 
provided only for the case when the diquark is assumed to be a composite 
object. 

We have to give some preliminary remarks before presenting the result. 
The Bialas-Bzdak model shows a singular behaviour at the dip position, 
having no real part in the amplitude. In order to give a meaning to the fit 
in this region 3 data points were left out from this dip region during the fit; 
these points are shown in red on the plots. 

Another important remark is that the TOTEM data covers the \ t\ range 
from 0.36 GeV up to 2.5 GeV and this range is applied in our minimization 
procedure to allow a comparison between the ISR and TOTEM results. 
Note that Bialas and Bzdak adjusted the model to account only for the 
slope and the value at t = together with the position of the minimum [T]. 
A different strategy is followed here, since the theoretical curve was fitted 
directly to the experimental data points using the CERN MINUIT package. 
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p+p p+p, diquark as a single entity,Vs=23.5 GeV 
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Fig. 11. Results of MINUIT fits at ISR energies when the diquark is assumed to be 
a single entity. The confidence level is higher than 0.1%, which means that the fit 
quality is acceptable. The fit was mode on the 0.36 - 2.5 GeV \t\ range according to 
[16[ . As the model is singular around the dip, 3 data points at the dip were left out 
from the fit. The parameter values are given with statistical errors. The systematic 
errors are not yet investigated and the effects of the correlation between model 
parameters are not yet determined, in this sense this result is still preliminary. 
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p+p p+p, diquark as a single entity,Vs=30.7 GeV 
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12. Same as Fig. [TTl but for the energy ^/s = 30.7 GeV. 
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p+p p+p, diquark as a single entity,Vs=52.8 GeV 
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13. Same as Fig. [TTl but for the energy ^/s = 52.8 GeV. 
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p+p p+p, diquark as a single entity,Vs=62.5 GeV 
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14. Same as Fig. [TTl but for the energy ^/s = 62.5 GeV. 
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p+p p+p, diquark as a single entity,Vs=7000.0 GeV 
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Fig. 15. The result of the fit at LHC at 7 TeV when the diquark is assumed to be 
one entity. The obtained parameter values are given with statistical errors. Note 
that CL is below 0.1%, so the quality of this fit is not acceptable. The systematic 
errors are not yet investigated and the effects of the correlation between model 
parameters are not yet determined, in these sense this result is still preliminary. 



4-2. The diquark is assumed to be a composite object 

In this section the MINUIT results are presented for the ISR [14^ [T5] 
and TOTEM [16j proton-proton elastic scattering data using the as- 
sumption that the diquark is a composite qq object. The results are illus- 
trated on Fig. [TBll2Ui The confidence levels, and model parameters with 
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0.26 ± 0.01 


0.26 ± 0.01 


0.30 ± 0.01 


0.41 ± 0.01 


Rd [fm] 


0.71 ± 0.01 


0.72 ± 0.01 


0.73 ± 0.01 


0.77 ± 0.01 


0.80 ± 0.01 


A 


0.68 ± 0.03 


0.55 ± 0.02 


0.56 ± 0.03 


0.60 ± 0.04 


0.38 ± 0.01 


^qq 


0.57 ± 0.01 


1.00 ± 0.97 


1.00 ± 0.80 


1.0 ± 0.66 


0.94 ± 0.04 


X'/NDF 


60.6/44 


21.5/32 


50.2/33 


43.73/33 


4355.9/75 


CL [%] 


4.9 


92.0 


2.8 


10.0 


0.0 


^elastic [t^kI 

^totai y"^"\ 


6.0 


5.6 


5.9 


8.7 


20.3 



Table 1. The overall fit quality and resulting parameters of the fit at ISR energies 
including the LHC result at 7 TeV. The diquark is assumed to be a single entity. 
The obtained parameter values are given with statistical errors. The systematic 
errors are not yet investigated and the effects of the correlation between model 
parameters are not yet determined, so these results in this sense are preliminary. 



their errors are summarized in Table El 



25 



p+p p+p, diquark as a qq entity,\fs=23.5 GeV 



3 points left out 
Fit range = 0.36 < t < 2.50 GeW^ 
X^/NDF = 64.42/44= 1.46 
CL = 2.397 % 



Data points 

Not used 
Fitted theory 
Extrapolation 




Fig. 16. The result of the fit at the 23.5 GeV ISR energy when the diquark 
is assumed to be a composite entity. The confidence level is higher than 0.1%, 
which means that the fit quality is acceptable. The fit was mode on the 0.36 - 2.5 
GeV \t\ range according to [16^. As the model is singular around the dip, 3 data 
points at the dip were left out from the fit. The parameter values are given with 
statistical errors. The systematic errors are not yet investigated and the effects 
of the correlation between model parameters are not yet determined, in this sense 
this fit result is still preliminary. 
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p+p p+p, diquark as a qq entity,\fs=30.7 GeV 



3 points left out 
Fit range = 0.36 < t < 2.50 GeV^ 

I NDF = 30.42/32= 0.95 
CL = 54.644 % 



o Data points 

• Not used 

Fitted theory 

Extrapolation 




-t [GeV^] 

Fig. 17. Same as Fig [1^1 except that it is provided for ^/s = 30.7 GeV. 



p+p p+p, diquark as a qq entity,\fs=52.8 GeV 
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3 points left out 
Fit range = 0.36 < t < 2.50 GeV^ 
X^/NDF = 64.09/33= 1.94 
CL = 0.094 % 



o Data points 

• Not used 

Fitted theory 

Extrapolation 



>i = 1.00 ±0.22 
A„„ = 0.63 ±0.03 



: 0.1 6 ±0.02 fm 
: 0.35 ±0.01 fm 



: 0.75 ±0.01 fm 



0.5 



1 



,1-5 



2.5 



-t [GeV^] 

Fig. 18. Same as FigUl given for = 52.8 GeV. 
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p+p p+p, diquark as a qq entity,\fs=62.5 GeV 



I 10^ 
E 1 



,-1 
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§10 



£10-' 

10^ 
10^ 
10-^ 
10^ 



3 points left out 
Fit range = 0.36 < t < 2.50 GeV' 
X^/NDF = 46.39/33= 1.41 
CL = 6.106% 




o Data points 

• Not used 

Fitted theory 

Extrapolation 



X = 0.01 ± 0.75 
A,q = 0.94 ±0.08 

Rp = 0.17±0.02 fm 
Rq = 0.32 ±0.01 fm 
Rh = 0.77 ±0.01 fm 




-t [GeV^] 

Fig. 19. Same as Fig [HI calculated for y/s = 62.5 GeV. 
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p+p p+p, diquark as a qq entity,\fs =7000.0 GeV 



I 10 
E 1 



,-1 



(0 

§10 



£10-' 
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10^ 



10' 



10" 
10' 



Preliminary resuit 
Fit range = 0.36 < t < 2.50 GeV^ 

I NDF = 4242.89/75= 56.57 
CL = 0.000 % 



o Data points 

Fitted theory 

Extrapolation 




X = 0.01 ± 0.01 
A„„ = 0.76 ±0.01 



: 0.54 ±0.01 fm 
: 0.43 ±0.01 fm 



: 0.73 ±0.01 fm 




t [GeV^] 



Fig. 20. Same as Fig [121 but for = 7 TeV. Note that this fit is not acceptable 
since CL is below 0.1%. 



^.3. Total cross sections in the composite diquark model 

The total inelastic cross sections for the quark-quark, quark-diquark and 
diquark-diquark subcollisions were analysed according to formula dH). The 
detailed results are collected in Table [3l while the average ratios for the 
described ISR energies are 



Tqq ■ Tqd ■(Tdd = '^- (1-93 ± 0.03) ! (3.64 ± 0.1), 



(15) 
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[GeV] 


23.5 


30.6 


52.9 


62.5 


7000 


Rp [fm] 


0.23 ± 0.01 


0.18 ± 0.02 


0.16 ± 0.02 


0.17 ± 0.02 


0.54 ± 0.01 


Rq [fill] 


0.28 ± 0.01 


0.34 ± 0.01 


0.35 ± 0.01 


0.32 ± 0.01 


0.43 ± 0.01 


Rd [fm] 


0.72 ± 0.01 


0.73 ± 0.01 


0.75 ± 0.01 


0.77 ± 0.01 


0.73 ± 0.01 


A 


0.62 ± 0.03 


1.00 ± 0.64 


1.00 ± 0.22 


0.01 ± 0.75 


0.01 ± 0.01 


^qq 


1.00 ± 0.93 


0.65 ± 0.03 


0.63 ± 0.03 


0.94 ± 0.08 


0.76 ± 0.01 


X'/NDF 


64.42/44 


30.42/32 


64.1/33 


46.4/33 


4242.9/75 


CL [%] 


2.40 


54.6 


0.1 


6.1 


0.0 


^elastic [t^kI 

^totai y"^°\ 


6.5 


6.0 


6.4 


8.4 


13.1 



Table 2. The overall fit quality and resulting parameters of the fit at the ISR 
energies including the LHC result at 7 TeV. The diquark is a assumed to be a 
qq entity. The obtained parameter values are given with statistical errors. The 
systematic errors are not yet investigated and the effects of the correlation between 
model parameters are not yet determined, so these results are preliminary. 



V~s [GeV] 


23.5 


30.6 


52.9 


62.5 


7000 


Cfqd/'^qq 


1.92 ± 0.08 


1.93 ± 0.01 


1.93 ± 0.01 


1.92 ± 0.01 


1.87 ± 0.01 


<^dd/<^qq 


3.64 ± 0.34 


3.66 ± 0.03 


3.67 ± 0.03 


3.61 ± 0.05 


3.38 ± 0.02 



Table 3. The ratios of the total inelastic cross sections for the quark-quark, quark- 
diquark and diquark-diquark processes for the ISR and LHC energies using the 
composite diquark hypothesis. 



which is close to the ideal 1: 2 : 4 ratio, confirming the assumption of having 
two quarks inside the diquark, amended with some shadowing which is 3.5% 
and 9% respectively. At 7 TeV the ratios are slightly different from (|15|) 

1 : (1.87 ± 0.01) : (3.38 ± 0.02), (16) 

which shows that the shadowing is stronger, 6.5% and 16% percent respec- 
tively. 

The total elastic scattering cross sections were also determined, as given 
in Table [2j Note, that the errors on the total cross-section have not yet been 
obtained reliably at this point. The reason is that the errors of total cross- 
sections are very strongly correlated with the errors of Aqq. However the 
current fits cannot determine precisely the value of this parameter: its nearly 
100 % relative error indicates the approximate insensitivity of the results 
on this value. Hence one has to study the possibility of fixing this and also 
the parameters like A to reasonable values and to see if in such scenarios the 
quality of the results remains the same or not. These studies and the final 
values of the total elastic scattering cross-sections will be submitted for a 
separate publication. 
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5. Conclusion and outlook 

A systematic study of fit quality as well as the fit parameters under 
similar circumstances has been performed for the Bialas - Bzdak model in 
a wide energy range from ISR to LHC energies. The model gives a good 
description of the ISR data, which means that the CL is acceptable on the 
ISR energies if 3 data points at the dip were left out from the fit. The total 
proton-proton cross section ( "size" of the proton) clearly seems to grow with 
energy but the model fails at this energy domain as CL is not acceptable 
at 7 TeV. This preliminary result does not include systematics, and the 
correlation between the fit parameters are still under investigation. 

An important shortcoming of the quark-diquark model of protons is that 
it ignores the real part of the elastic scattering amplitude. This leads to a 
singular behaviour at the diffractive minimum, which is apparently a more 
and more serious model limitation with increasing the energy of the p+p 
collisions. 

The evaluated ratios of the quark-quark, quark-diquark and diquark- 
diquark total inelastic cross-sections were found to deviate more and more 
from the ideal 1:2:4 ratio with increasing energies. In the ISR energy 
range the deviations from this ideal value were not yet significant, indicating 
lack of significant shadowing effects. However at the current LHC energy 
of ^/s = 7 TeV, a significant decrease compared to these ideal ratios were 
found, which possibly may indicate an increased role of shadowing at CERN 
LHC energies. 

The final conclusions of this study will be summarized separately and 
submitted for a publication. The current status corresponds to the level of 
our understanding around September 2011, at the time of the 2011 Summer 
School on Diffraction in Heidelberg, where these results were first presented. 

Finally let us note that the TOTEM Collaboration extended recently 
the measurement of the differential elastic p+p scattering cross-sections to 
low values of |t| in ref. ^17j, allowing one to extrapolate to the optical point 
at t = and to determine the total elastic and the total scattering cross- 
sections of p+p collisions at ^/s = 7 TeV for the first time, but these data 
were yet not utilized in our analysis. 

As a final remark, at the time of the completion of this conference con- 
tribution it is inspiring to see the great theoretical interest in the differential 
elastic cross-section measurement of TOTEM at LHC, as evidenced by the 
increasing amount of theoretical interpretations and successfull descriptions 
of certain aspects of these data. As our contribution is not intended to be 
a review on the interpretation of TOTEM data, we just would like to call 
attention to some of the most interesting approaches of describing these 
measurements as evidenced in refs. [Ml [H [2U1 [JH [22l [23l [Ml [25] . 
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7. Appendix 

Two of the Dirac 6 functions in ([T|) induce the following transformation 



in the transverse diquark and quark position variables 

Sd = — ASq, Sd' = —XSq'. 



(17) 



Hence four Gaussian integration remain, which lead us to the following 
result 



X e 



^V^ _u2 r 
— e n 



where the coefficients Cab are abbrevations, and 

VL= 4^; + (1 + A)^ {Cqd + Cdq) [v + Cqq + A^Crfrf] + 



while 



r 



+ 



+ (1 - A)^ V {Cqd + Cdq) + (1 + A)^ CqdCdq 



4U + (1 + A)^ {Cqd + Cdq) V {Cqq + Cdd) + (1 + A)^ CqqCdd 
4u + (1 + A)^ {Cqq + Cdd) V {Cqd + Cdq) + (1 + A)^ CqdCdq 



+ 



(18) 



(19) 



(20) 
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